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Abstract 

The magic proton and neutron numbers are searched in the superheavy region with proton 
number Z=100 - 140 and neutron number N= (Z+30) - (2Z+32) by the relativistic continuum 

m 

Hartree-Bogoliubov (RCHB) theory with interactions NL1, NL3, NLSH, TM1, TW99, DD-ME1, 



PK1, and PK1R. Based on the two-nucleon separation energies S^p and S^m the two-nucleon gaps 
&2p and 62m the shell correction energies E p shell and E™ hell , the pairing energies E^ air and E£ air , and 
the pairing gaps A p and A n , Z=120, 132, and 138 and iV=172, 184, 198, 228, 238, and 258 are 
suggested to be the magic numbers within the present approach. The a-decay half-lives are also 
discussed. In addition, the potential energy surfaces of possible doubly magic nuclei are obtained by 
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the deformation-constrained relativistic mean field (RMF) theory, and the shell effects stabilizing 
the nuclei are investigated. Furthermore, the formation cross sections of f 72 120 and 1^120 at the 
optimal excitation energy are estimated by a phenomenological cold fusion reactions model with 
the structure information extracted from the constrained RMF calculation. 

PACS numbers: 21.10.Dr, 21.60.Jz, 24.10.Jv, 24.60.Dr, 27.90.+b 
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I. INTRODUCTION 



The structure of heavy and superheavy nuclei has been an interesting field of nuclear 
physics research during the last decades. Since the stability of superheavy nuclei is mainly 
determined by the shell effects, it is important to find out the regions in the (Z, N) plane 
where the shell effects are strong and the long lifetime of the superheavy nuclei can be 
expected. Exploring the limit of nuclear charge and mass is a long-term goal of nuclear 
physics, and searching for the superheavy nuclei has been enlivened by the hope of creating 
nuclei with masses and charges much larger than those we are familiar with. The limit on 
stability in superheavy nuclear region is essential for understanding not only the nuclear 
structure, but also the structure of the stars and the evolution of the universe. In the 
classical droplet theory, the superheavy nuclei can not exist due to their large Coulomb 
potentials. It is the shell effects that play a major role for the very existence of nuclei with 
magic numbers and provide higher stability and also higher abundance as compared to their 
neighbors. A semblance of the same would work also for superheavy nuclei, if there were 
magic numbers in this region. Consequently, these nuclei will be guarded against a faster 
decay by fission as compared to their non-magic counterparts, and have more opportunities 
to be bound. 

The pioneer theoretical work on the superheavy elements can be found in the 1960's 
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]. A series of calculations based on the macroscopic-microscopic method (Nilsson- 
Strutinsky approach) with the folded- Yukawa deformed single-particle potential [5] and with 
the Woods-Saxon deformed single-particle potential ((J UQ| are successful in reproducing 
the well-known a-decay half-lives of heavy elements. The existence of island of superheavy 
nuclei would be separated in neutron and proton number from the known heavy elements by 
a region of much higher instability. Experimentally, the superheavy element s up to Z=116 
are synthesized or claimed to be synthesized Q, H Q 13, 14, 15, H H fl Q by 



different heavy ion reaction types including the cold fusion reactions. The usual view is 
that although these isotopes are very heavy indeed, they are not examples of the originally 
sought island of superheavy elements. 

There is no consensus among different theories with regard to the center of the island 
of superheavy nuclei. Based upon phenomenological models such as finite-range droplet 
model (FRDM), the shell closure were predicted at Z—114 and iV=184 j2Q|. Additionally, 



2 



he FRDM also predicts larger shell gaps at Z=10A, 106, 108, 110 and at iV=162, 164 
2o| . Based on Nilsson-Strutinsky scheme, a similar pattern of deformed nuclei have been 



predicted about Z=108 and iV=162 as in FRDM [8 

the question whether the macroscopic approaches which apply to the region of /5-stabi 
line can be extrapolated to the superheavy nuclei. Recently, microscopic calculations 



2l|. However, the main obstacle is 
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28, 



22, 



291 ] are also attempted in describing the superheavy nuclei. In the 



framework of relativistic Hartree-Bogoliubov (RHB) theory, calculation with a finite-range 
pairing force of Gogny interaction Dl and effective interaction NLSH show that Z=114 and 
iV=160, iV=166, iV=184 exhibit stability compared to their neighbors and indications for a 
doubly magic character at Z=106 and iV=160 are also observed 23]. The Skyrme Hartree- 
Fock (SHF) method with interactions SkP and SLy7 predicts magic numbers at Z=126 and 
iV=184, and also predicts the increased stability due to the deformed shell effects at iV=162 



24j. Considering non-relativistic SHF effective interactions SkM*, SkP, SLy6, Skll, SkI3, 
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4 and relativistic mean field (RMF) effective interactions PL-40, NLSH, NL-Z, TM1, Ref. 
J gave the doubly magic spherical nuclei 184114, 172I2O and 184126 based on two-nucleon 
gaps 52 P and 52 n - The common point of such studies is that the magic numbers depend on 
effective interactions. The study of superheavy nuclei remains a challenge for interaction, 
which have to be rigorously tested by a wide variety of nuclear properties throughout the 
periodic table. 

The extrapolation towards superheavy nuclei challenges the predictive power of nuclear 
structure models. The relativistic mean field (RMF) theory has achieved a great success in 

and it is believed as one of the best candidates for 
application to superheavy nuclei. The relativistic continuum Hartree-Bogoliubov (RCHB) 
theory is the extension of the RMF and the Bogoliubov transformation in the coordinate 
representation As the RCHB formalism allows for the proper description of the coupling 
between the bound states and the continuum by the pairing force, therefore it is suitable 
not only for stable nuclei but also for the nuclei near the drip line 

The 

RCHB theory has shown a remarkable success in the description of nuclei with unusual N/ Z 
ratio, e.g., the halo in 11 Li|32] and the description of the giant halo at the neutron drip line 
in Zr and Ca isotopes 33J, |36j . A systematical study on 259 Db and its a-daughter nuclei were 
also achieved J^. In principle, only a calculation in a large multidimensional deformation 
space can definitively decide the appropriate ground-state shape. However, as the traditional 
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superheavy nuclei are expected to be spherical and are located on the nuclear chart around 
a spherical doubly magic nucleus next to iPJPb (37], the RCHB theory with the assumption 
of spherical shape can be applied to this preliminary scan for magic numbers. It is the aim 
of this paper to predict doubly magic superheavy nuclei within the RCHB theory using a 
variant of the RMF interactions. As most of the presently found superheavy nuclei have 
large deformations, we will study the potential energy surface with the deformed-constrained 
RMF theory in order to confirm that the spherical configuration exists in these doubly magic 
nuclei. 

The material contained in this study is organized as follows. The framework of RCHB 
and constrained RMF theory are outlined in Section II. Section III contains the predictions 
made from the two-nucleon separation energies S 2p and S 2n , the two-nucleon gaps 5 2p and 
6 2n , the shell correction energies E p shell and E™ hell , the pairing energies E v pair and E^ a ; r , and 



the pairing 



;ap A p and A n , obtained in the RC 



TB theory with interactions NL1 |38[, NL3 



NLSH 4Q|, TM1 Q, TW99 42 1, DD-ME1 3> PK1 and PK1R Q The a-decay half- 
lives are also discussed there. The binding energy and the corresponding macroscopic energy 
as a function of deformation obtained in the deformation-constrained RMF calculation for 
.2=120 are presented in Sec. IV. Additionally, the formation cross sections for the possible 
doubly magic nuclei 292 1 20 and 304 120 at the optimal excitation energy are estimated by a 
phenomenological model proposed by Smolahczuk j3] with the input from the constrained 
RMF theory. Finally, Section V contains the main conclusions of this work. 

II. FRAMEWORK 

A. Relativistic continuum Hartree-Bogoliubov theory 

The basic ansatz of the RMF theory is a Lagrangian density where nucleons are described 
as Dirac particles which interact via the exchange of various mesons. The Lagrangian density 
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considered has the form: 



L = if} 
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-\f» v f» v (1) 

where ip is Dirac spinor and ip = V f 7°- M , 

m ff , m^, and m p are the nucleon, a, to and 
p meson masses respectively, while g a , g 2 , #3, g u , C3, g p , d 3 , and e 2 /47r = 1/137 are the 
corresponding coupling constants for the mesons and the photon. The field tensors of the 
vector mesons (u and p mesons) and of the electromagnetic fields take the following form: 

W = d"u v - d u u» 

W v = - SP? - 2g p ? x ? (2) 
= d»A u - d u A» 

With the classical variational principle, one can obtain the coupled equations of motion, 
which contain the Dirac equation for the nucleons and the Klein-Gordon type equations for 
the mesons and the photon. The coupled equations are self-consistently solved by iterations. 

For the RMF theory, however, as the classical meson fields are used, the equations of 
motion for nucleons derived from Eq. ((TJ do not contain pairing interaction. In order to 
have two-body interaction, one has to quantize the meson fields which leads to a Hamiltonian 
with two-body interaction. Following the standard procedure of Bogoliubov transformation, 
a Dirac Hartree-Bogoliubov equation could be derived and then a unified description of the 
mean field and pairing correlation in nuclei could be achieved j^. The Dirac Hartree- 
Bogoliubov equations are as following: 

= eM, (3) 




where, h = a ■ p + g u u + g p r 3 p + (3(M + g a a) is the Dirac Hamiltonian and the Fock term has 
been neglected as is usually done in RMF. A Lagrange multiplier A is introduced to adjust 
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the proper particle number of neutrons and protons. The pairing potential A in Eq. (J3J) is : 



£ ( 4 ) 



A ab 2 

cd 

It is obtained from one-meson exchange interaction ^6cd m the pp-channel and the pairing 
tensor k = ipljipy- The interaction V a b c d used for the pairing potential in Eq. (J1J is either the 
density-dependent two-body force of zero range or Gogny-type finite range force as shown 
in Ref. 3jJ. Considering the consuming computational time, we prefer to use the density- 
dependent zero range 5 force with the interaction strength Vq and the nuclear matter density 
Po as 

V(n, r 2 ) = V 5(n - f 2 )\ [1 - a x ■ a 2 ] (l - ^) . (5) 

4 V Po / 

The relativistic continuum Hartree-Bogoliubov theory solves the RHB equations in the 
coordinate representation in order to describe both the continuum and the bound states self- 
consistently 31]. It is then applicable to not only stable nuclei but also exotic nuclei. When 
the 5-force in the pairing channel is used, the RCHB equations are a set of four coupled 
differential equations, which are solved in a self-consistent way by t he shooting method and 
the Runge-Kutta algorithm. More details can be found in Ref. |3jJ. 

In this paper, the RCHB equations are solved in a spherical box with radius R=25 fm, 
the step size of 0.1 fm, and proper boundary conditions. It has been shown in Refs. 0, 
that the results do not depend on the box size when R > 15 fm. The strength Vq of the 
pairing force for proton and neutron is fixed at -650 MeV-fm -3 . For p we use the nuclear 
matter density 0.152 fm -3 . 



B. Deformation-constrained RMF theory 

To examine the ground state geometrical configuration of certain heavy or superheavy 
nucleus, as well as its fission barrier, one has to obtain the potential energy surface. The 
potential energy curve can be calculated microscopically by the deformation-constrained 
RMF theory. The binding energy at certain deformation value is obtained by constraining 
;he quadrupole moment {Q 2 ) to a given value [i 2 in the expectation value of the Hamiltonian 

{H') = (H) + ^C,((Q 2 )-fi 2 ) 2 , (6) 
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where is the constraint multiplier. The expected value of the quadrupole moment oper- 
ator is < Q2 >=< Q2 > + < Q2 >> where the quadrupole moment operator for neutrons 
(protons) is calculated by < Q% >=< 2r 2 P 2 (cos6) > n ( P )- The equations of motion are 
solved by expanding the Dirac spinors in harmonic oscillator basis in the cylindrical coor- 
dinates. The size of the basis is crucial and responsible for the reliability of the result. By 
varying /ig, the binding energy at different deformation can be obtained. The pairing is 

taken into account by the constant gap approximation (BCS) in which the pairing gap is 
12 

taken as —= for even number nucleons. An enormous advantage of microscopic constrained 
VA 

RMF calculations compared with the phenomenological calculations is that in the micro- 



scopic constrained calculation we need only follow a one-parameter line, e.g. fa, whi 



the phenomenological methods we must calculate multidimensional energy surfaces 47 1 
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III. MAGIC NUMBERS IN SUPERHEAVY NUCLEI 



The doubly magic nuclei are searched for the even-even nuclei (more than 1200 nuclei) 
in the nuclear chart with proton number Z=100 - 140 and neutron number N= (Z+30) - 
2Z+32) (i.e. , iV=130 - 312) by the RCHB calculation with the effective interactions NL1 



PK1 and PK1R U4. These 



38j, NL3 [39|, NLSH 0, TMl |41|, TW99 |42j, DD-ME1 
nuclei are presented as shaded area in Fig. The long-dashed line and solid line in Fig. ^ 
represent the /3-stability line of the parametrization 



N - Z = 6 x 10~ 3 • A 5/3 



from Ref. 



. 3 



and 



(7) 



(8) 



from Ref. 



1.98 + 0.0155 • A 2 / 3 

49J, respectively. The black dots lying at the lower-left corner represent the 
superheavy nuclei which have been observed or declared to be observed experimentally j3 • 
The vertical and horizontal lines in Fig. ^represent the possible magic proton and neutron 
numbers which will be discussed later. 

For RMF calculations, there exist many effective interactions, which provide nearly equal 
quality of description for stable nuclei. However, when they are applied to the exotic nuclei, 
the differences appear. In view of the uncertainties, the effective interactions considered here 
are four frequently used interactions NL1, NL3, NLSH, TMl, density-dependent interactions 



TW99 and DD-ME1, as well as the newly developed PK1 and PK1R which were adjusted 
with the elaborate consideration of the cent-of-mass correction. The effective interactions 
PK1 and PK1R are listed in Table □ 

In the following four subsections, the signatures of shell closure on the two-nucleon sep- 
aration energies S 2p and S 2n , the two-nucleon gaps S 2p and S 2n , the shell correction energies 
E shell an d E shell? the pairing energies E v pair and E pair and the pairing gaps A p and A n ob- 
tained in RCHB calculation will be analyzed consecutively. The a-decay half-lives T a will 
be discussed afterwards. At last, the prediction of the magic proton and neutron numbers 
in superheavy nuclei will be summarized. 

A. Two-nucleon separation energies 

The two-nucleon separation energy, 

S 2n (N, Z) = Eb(N, Z) -E b (N- 2, Z), S 2p (N, Z) = E B (N, Z) - Eb(N, Z-2) (9) 

is better to quantify shell effects than the single-nucleon separation energy due to the absence 
of odd-even effects, here Eb(N, Z) is the binding energy of the nucleus N+Z Z. For an isotonic 
chain, the S 2p become smaller with the proton number Z. A jump of S 2p indicates the 
occurrence of a proton shell closure. 

For the sake of a clear presentation of the results, the two-proton separation energies S 2p 
from the RCHB calculation for even-even nuclei with .2 = 102 - 140 as a function of mass 
number A with the effective interactions NL1, NL3, NLSH, TM1, TW99, DD-ME1, PK1, 
and PK1R are shown in Fig. |21 In the figure, S 2p (N, Z) for the same Z are connected to 
form a curve so that the jumps of S 2p compose a gap between the curves. The 20 S 2p curves 
in each subfigure correspond to isotopic chains for Z ranging from 102 to 140 with a step of 
two. The curve at the upper- left corner is for Z = 102 and the neighboring one is for Z=104, 
and so on. 

For an isotonic chain, the S 2p gradually decrease from about 25 MeV to zero for the 
curve at the upper-left corner to the lower-right corner. Each curve is roughly in parallel 
with the other. A gap is enchased if a magic proton number exists. Take S 2p for the effective 
interaction NL1 as an example: the gap between Z=120 and Z=122 (abbr. at Z=120) is 
large and stable, the gap at Z=132 is visible just in the region iV=220 - 250, and the gap 
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at Z=138 is relatively large but unstable. The gap at Z=138 becomes weak at 240. 

For all effective interactions, the magic proton numbers Z=120, 132, and 138 are common 
while Z=106 is observed only for NL3, NLSH, TW99, PK1, and PK1R. Generally, the gaps 
do not exist for all isotopes, e.g., the gap at Z=132 for NL1 exists only for iV=220 - 250, 
and the gap at Z=138 for PK1 exists only for iV=230 - 270. 

Fig. El is the same as Fig. |21 but for S^n with iV=132 - 312. Each curve corresponds 
to an isotonic chain. The upper- left curve is for iV=132 and the difference of the neutron 
number between the neighboring curves is 2. 

The gaps of S2 n are more complicated than those of S 2p - For fixed Z, S 2n decrease 
gradually with N. For effective interaction NL1, there are nine gaps at iV=138, 172, 184, 
198, 228, 238, 252, 258, and 274. The size and the shape of the gaps differ notablely: The 
gaps at iV=172, 184, and 258 are relatively larger than the other gaps. 

For all effective interactions, S 2n show distinguishable gaps at iV = 138, 172, 184, 198, 
228, 238, 252, 258, and 274. Such gaps mark the magic neutron number. Apart from 
the above mentioned common gaps, there are also other gaps which appear only for some 
effective interactions, i.e., N=164 for NLSH, TW99, DD-ME1, PK1 and PK1R, and N=216 
for NLSH, TW99, PK1, and PK1R. For the magic numbers marked either by the common 
gaps or the interaction-dependent gaps, the shell quenching phenomena, i.e., the gaps at 
JV=184, 198 (NL1), 216(NLSH, TW99, PK1, and PK1R), 228 (NL1, NL3, NLSH, TM1, 
PK1, and PK1R), 238, 252, and 258 (NLSH, TM1, PK1, and PK1R) appear only for certain 
Z, are observed. 

Summarizing the above results, we may say that based on two-nucleon separation energies 
S 2p and S 2n , Z=120, 132, and 138 ,and JV=138, 172, 184, 198, 228, 238, 252, 258, and 274, 
are the magic numbers predicted in the RCHB calculation. They are independent of the 
effective interactions used in the calculation. There are also some other magic numbers, e.g., 
Z— 106 and iV= 164 and 216, which appear only for some effective interactions. 
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B. Two-nucleon gaps 



The changes of the two-nucleon separation energies can also be quantified by the second 
difference of the binding energies, i.e., the two-nucleon gaps: 

5 2n (N, Z) = 2E B (N 1 Z) - E B (N + 2,Z)- E B (N -2,Z) = S 2n (N, Z) - S 2n (N + 2, Z}0) 
5 2p (N, Z) = 2E B (N, Z) - E B (N, Z + 2) - E B (N, Z - 2) = S 2p (N, Z) - S 2p (N, Z + tfLl) 

A peak of the two-nucleon gap implicates the drastic changes of the two-nucleon separation 
energies, which can be used as the evidence of the magic number existence. 

The two-proton gaps 5 2p from the RCHB calculation for even-even nuclei with Z=102 
- 138 as a function of Z with the above eight effective interactions are shown in Fig. 0] 
The two-proton gaps 5 2p with the same N are connected as a curve. A peak at certain Z 
in the curve suggests the existence of magic proton number. The sharpness of the peaks 
represent the goodness of the magic numbers while the quenching effects are associated with 
the bundle of the curves at the certain Z. The size of the gaps of two-proton separation 
energies S 2p in Fig. |2] correspond to the magnitude of the peaks of two-proton gaps S 2p 
in Fig. E] . Generally the magic numbers indicated by two-proton separation energies S 2p 
are supported by two-proton gaps 5 2p . It is observed that common magic proton numbers 
Z=120, 132, and 138 exist for all effective interactions while Z=106 is observed only for NL3, 
NLSH, TW99, DD-ME1, PK1, and PK1R. Furthermore, the peak at Z=1U for NLSH and 
TW99 and the peak at Z=12Q for NL1 are also observed, though they are not so obvious as 
that at Z=120. 

As counterpart of two-proton gaps S 2p , two- neutron gaps 5 2n from the RCHB calculation 
for even-even nuclei iV=132 - 308 with effective interactions NL1, NL3, NLSH, TM1, TW99, 
DD-ME1, PK1, and PK1R are shown in Fig. There are interaction- independent peaks of 
two-neutron gaps 5 2n at iV=138, 172, 184, 198, 228, 238, 258, and 274 which are consistent 
with the common gaps of two-neutron separation energies S 2n in Fig. El Moreover, small 
peaks are also observed for iV=154 (NLSH and TW99), iV=164 (NL3, NLSH, TW99, DD- 
ME1, PK1, and PK1R), and iV=216 (NLSH, TW99, PK1, and PK1R). 

In general, based on two-nucleon gaps 5 2p and 5 2n , the magic numbers Z=120, 132, 
and 138, and JV=138, 172, 184, 198, 228, 238, 258, and 274 can also be observed for RCHB 
calculation with all effective interactions. Apart from Z=106 and iV=164 and 216, additional 



10 



magic numbers Z=114 and 126, iV=154 appear in the two-nucleon gaps 52 P and #2n in Figs. 
Eland El for some effective interactions. 



C. Shell correction energies 

The shell correction energy, representing the overall behavior of the single-particle spectra, 
may be a good candidate to identify the magicity. This quantity, which is derived from the 
single-particle spectra, is defined as the difference between the total single-particle energy 
E and the smooth single-particle energy E: 

N(Z) a 

E shell = E-E=J2 e i~ 2 / eg(e)de. (12) 

i=1 -oo 

where N(Z) is the particle number, is the single-particle energy, A is the smoothed Fermi 
level and g(e) is the smoothed level density. The shell correction energy provides an indicator 
about the deviation in the level structure of nuclei away from uniformly distributed ones. 
A large negative shell correction energy corresponds to shell closure at particular nucleon 
number. 

The shell correction energy can be obtained by the Strutinsky procedure, in which the 

A 

smoothed Fermi level A is determined by the particle number equation N(Z) = 2 J g(e)de. 

— oo 

The smoothed level density g(e) takes the form: 

OO / \ 

9(e) = - / HP 5(e' - e 4 ) /(— )cfe' = ± £ /(^). (13) 

1 _ 2 

where 7 is the smoothing range, the folding function f(x) = —^e x P(x), and P(x) is an 

V 71 " 

associated Laguerre polynomial L l J 2 (x 2 ). 

In present calculation, the smoothing range is 7 = 1.2 and the order of the associated 

Laguerre polynomial is s = 3. The unit of 7 is 41 A^ 1 ^ 3 ^1 ± - — — — j MeV where the 

plus (minus) sign holds for neutrons (protons). Then the particle number equation can be 
evaluated as: 

N{Z) = 2 f; - ^=(57 tl - m\ + At!)}, (14) 

i=i * 
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where ti = — and erf(tj) is the error function. The corresponding smooth single-particle 

7 

energy E takes the form: 
E = 2 f; {e ^-f^ - ^(57U - 32tf + 4tf)] - 7^4(15 - 90^ 2 + 60*J - 8tf)}. 

i=l 

(15) 

With the proton single-particle spectra in canonical basis from the RCHB calculation for 
even-even nuclei with Z=100 - 140 as a function of Z with effective interactions NL1, NL3, 
NLSH, TM1, TW99, DD-ME1, PK1, and PK1R, the shell correction energies for proton 
E p shell are shown in Fig. El The nuclei with the same N are connected as a curve. A deep 
valley at certain Z on a curve hints the magicity. The valleys at Z=120, 132, and 138 are 
conspicuous for all effective interactions. Possible shell closure at Z=106 for NL1, NL3, 
NLSH, TW99, DD-ME1, PK1, and PK1R, at Z=1U for NLSH, TW99, PK1, and PK1R, 
and at Z=12Q for NL1 are also observed. The shell closure from the shell correction energies 
for proton E p shell in Fig. El are consistent with the peaks of the two-proton gaps <$2 P in Fig. 
El though the magnitude of magicity slightly differs. Similar to two-proton gaps <$2 P , the 
spread of the valleys of the shell correction energies for proton E p shell indicates the quenching 
of the magic number. In view of the shell correction energies for proton E P hell , the quenching 
phenomena of magic proton numbers in RCHB calculation are universal. These facts suggest 
that the magic proton numbers depend also on neutron number N. 

The corresponding shell correction energies for neutron E™ hell from the RCHB calculation 
for even-even nuclei with iV=130 - 312 as a function of N are demonstrated in Fig. [7| 
As the neutron numbers extends from iV=130 to iV=290, the valleys for shell correction 
energies for neutron E™ hell in Fig. are not so obvious as those for shell correction energies 
for proton E p shell in Fig. El However, similar conclusions as two-neutron separation energies 
Si n in Fig. El and two-neutron gaps bin in Fig- El for magic neutron numbers can also be 
drawn. For all effective interactions, iV=184 and 258 are found to be the common magic 
„» mb e rS . Slml .a r eolation has a, 8 o be done in Re, Q with SHF and RMF. PVo m the 
shell correction energies for neutron, although the magic neutron numbers are not so sharp 
as those got from two-neutron separation energies Sin and two- neutron gaps S^n, however, 
the small irregularity in the shell correction energies still correspond to the magic number 
discussed above. 

The first valley locates in iV=184 for all effective interactions and involves with iV=172 
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(NL3, NLSH, TM1, TW99, DD-ME1, PK1, and PK1R), even with iV=164 (NLSH and 
TW99). And the other valley is the mixture of the valleys at iV=228 (TW99 and DD-ME1), 
iV=252 (NLSH, TM1, PK1, and PK1R) and iV=258. These fine structures suggest the 
smearing of magic neutron numbers. 

Summarizing the above results, we may say that based on shell correction energies, the 
shell closures are smeared than those for two-nucleon separation energies S^ p and S2 n and 
two-nucleon gaps $2 P and 52 n , however we can observe magic numbers at Z=120, 132, and 
138 are common and strongly quenched, the shell closures near iV=172 and iV=184, and 
iV=228, iV=252, and iV=258 are blurred. The other magic numbers from the shell correction 
energies, such as 2 = 106, 114, and 126, appear only for some effective interactions. 



D. Pairing energies and pairing gaps 

In this subsection we will investigate the pairing effects, which can also provide reliable 
information of nuclear magic number, including the pairing energy and the effective pairing 
gap. The pairing energy in RHB theory is given by 

E pair = —TtAk, (16) 

where A and k are the pairing potential and the pairing tensor respectively. Generally 
speaking the pairing energies vanish at the closed shell and have a maximum values in the 
middle of two closed shells. Similarly, the effective pairing gap defined as 

£A^ 2 (2j + l) 
An(p) = N(z) , (17) 

can also be used to denote the shell closure, where the sum i is over the single particle 
,eve, S of _ or proton ln the eanonieal basifi 2j + ! is the de generacy number of 
the corresponding energy level i and N(Z) is the particle number, Aj is the pairing gap for 
level % in canonical basis determined by (E^ — A) ■ luiVi = Aj(t>? — uf) where A is the Fermi 
energy and the particle and hole amplitude u iy Vi satisfies uf + vf = 1. 

The pairing energies for proton E v vair and the effective pairing gaps for proton A p from the 
RCHB calculation for even-even nuclei with Z=100 - 140 as a function of Z with effective 
interactions NL1, NL3, NLSH, TM1, TW99, DD-ME1, PK1, and PK1R are shown in Figs. 
|H]and EH respectively. Both figures has the same pattern while the shell closures in Fig. E]are 
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more clear. It is observed that common magic proton numbers Z=120, 132, and 138 exist 
for all effective interactions while Z=106 exists only for NL1, NL3, NLSH, TW99, DD-ME1, 
PK1, and PK1R, Z=1U exists only for NL3, NLSH, TW99, DD-ME1, PK1, and PK1R, and 
Z=126 exists only for NL1, TW99, DD-ME1. Generally the shell closures are quenched for 
all magic proton numbers. It should be noted that the peak at Z=120 for NL1, TW99, and 
DD-ME1 is stable for all isotopes while for NLSH, TM1, PK1, and PK1R, it is quenched. 

Fig. El shows the pairing energies for neutron E™ air from the RCHB calculation for 
even-even nuclei with A=130 - 312 as a function of N with eight effective interactions. The 
vanishing of pairing at iV=138, 164, 172, 184, 198, 228, 238, 252, 258, and 274 are common 
for all interactions while the vanishing of pairing at A=154 and iV=216 exist only for NL3, 
NLSH, TW99, PK1, and PK1R, and while the vanishing of pairing at A=210 exists only 
for DD-ME1. The same pattern is observed by the effective pairing gaps for neutron A n in 
Fig. EU 

From the above discussion, the magic numbers Z=120, 132, and 138 and A=138, 164, 
172, 184, 198, 228, 238, 252, 258, and 274 are supported by the pairing energies and the 
effective pairing gaps. 

E. Alpha decay half-lives 

Alpha decay is one of the most predominant decay modes for superheavy nuclei. For an 
area of enhanced stability, the a-decay half-lives are expected to be longer than its neighbors. 
The phenomenological formula of Viola and Seaborg |52j for calculation of a-decay half-lives 
with the parameter in Ref. 

log 10 T a = (1.66175Z P - 8.5166)Q^ 1/2 - (0.20228Z P + 33.9069) (18) 

are used, where Z P is the proton number of a parent nucleus, Q a is the a-decay energy 
Q a = E B (Z - 2, N - 2) - E B {Z, N) + 28.295673 in MeV, and T a is in seconds. It should 
be noted that Eq. (jl8)l is based on the WKB approximation, and provides only a rather 
crude estimate of T a since it disregards many structure effects such as deformation, and 
configuration changes, etc. 

In Fig. ^1 we plot the half-lives T a from the RCHB calculation in logarithm scale as a 
function of neutron number N with effective interactions NL1, NL3, NLSH, TM1, TW99, 
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DD-ME1, PK1 and PK1R. The half-lives correspond to Ins, 1/is, 1ms, Is, lh, ly, and Iky 
are marked by the dashed lines. Each curve in this figure corresponds to an isotopic chain 
in the region 2=102-140. In Fig. the half-lives logioT a increase with N. The jumps of 
the curves correspond to magic proton numbers Z=106, 120, 132, and 138, which depend on 
the effective interactions. The peaks of each curve correspond to magic neutron numbers, 
i.e., A=172, 184, 228, 238 and 258. It can be seen that the magic numbers suggested by 
two- neutron separation energies 5*2n (two- neutron gaps #2n) or two-proton separation energies 
S2 P (two-proton gaps 5 2p ) can also be found here. 

So far, the RCHB theory with eight different interactions, namely NL1, NL3, NLSH, 
TM1, TW99, DD-ME1, PK1 and PK1R has been applied to the superheavy region with 
proton number Z=100 - 140 and neutron number N= (Z+30) - (2Z+32). After detailed 
analysises on two-nucleon separation energies and S2„, two-nucleon gaps b<i%> and b^m 
the shell correction energies E p sheU and E^ hell , as well as the pairing energies E v pair and E pair) 
the effective pairing gaps A p and A n , the magic numbers are predicted. Table ITD and Table 
IIHI list the possible magic proton and neutron numbers suggested by the above quantities 
and interactions. In these tables, the physical quantity (row-wise) corresponding to the 
magicity of certain proton number (column-wise) is marked as a black square for each 
interaction. According to Table HP and Table ITTT) the doubly magic superheavy nuclei can be 
the combination of magic proton numbers Z=120, 132, and 138 and magic neutron numbers 
iV=172, 184, 198, 228, 238, and 258. Such magic numbers are presented as gray lines in 
Fig. [TJ It should be noted that the crosses of lines in Fig. [T] do not necessary mean doubly 
magic nuclei as there are also shell quenching effects. However, these crosses include all the 
possible doubly magic nuclei in the state-of-the-art relativistic approach. 

IV. DOUBLY MAGIC NUCLEI IN SUPERHEAVY NUCLEI 

The magic numbers obtained before with the RCHB theory are based on the assumption 
of spherical geometrical configuration. However, the assumption is not always true for 
superheavy nuclei. In fact, most superheavy nuclei found experimentally are known to 
be deformed. It is worthy to investigate the potential energy surfaces in order to see the 
validity of spherical configuration. For this purpose, the deformation-constrained RMF 
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calculations have been done. The deformation parameter (3 2 of the harmonic oscillator basis 
(with 20 shells) is set to the expected deformation to obtain high accuracy and reduce the 
computing times. Although such calculation is very lengthy for superheavy nuclei, we have 
systematically calculated the potential energy surfaces for all the possible doubly magic 
nuclei with the effective interaction NL3. The calculation with other interactions for some 
typical doubly magic nuclei are also done for comparison. To save space, only the results 
for Z=120 isotopes with NL3 are presented here. 

In Fig. UJ the potential energy surfaces for nuclei 292 1 20 , 304 1 20 , 318 120 , 348 1 20 , 358 1 20, 
and 378 1 20 are presented with the solid lines in corresponding subfigures. These nuclei are 
clearly shown as the crosses with Z = 120 in Fig. ^ Although these nuclei have the same 
proton number, their potential energy surfaces are quite different from each other. For nuclei 
292 1 20 , 304 1 20 and 378 1 20, there is an obvious local minimum with the spherical configuration 
02 ~ 0, while another local minimum with large deformation (3 2 ~ 0.6 can be also clearly 
seen. For 318 120, the spherical local minimum is very shallow, and for 348 1 20 and 358 1 20, 
the spherical minimum is hardly seen and a local minimum with f3 2 ~ 0.25 appears instead. 
In addition to the local minima discussed, we should make remark on the absolute minima 
for these nuclei. For 292 1 20 and 3 120, the two minima with the different deformation 
are almost with the same energies, i.e. the so-called "shape coexistence", may exist. In 
particular, the spherical minimum is indeed the absolute minimum for 292 1 20 (Eb = -2064.3 
MeV with (3 2 ~ vs. E B = -2063.2 MeV with (3 2 ~ 0.6), while the ground state of 378 1 20 is 
quite delicate with E B = -2398.7 MeV for fa ~ and E B = -2398.9 MeV for f3 2 ~ 0.6. For 
304 1 20, the absolute minimum at (3 2 ~ 0.6 is much deeper (~ 6.1 MeV) than the spherical 
configuration. In addition for 348 1 20, the absolute minimum lies at (3 2 ~ 0.2, and it is a 
well-deformed nucleus according to the present calculation. 

To investigate the role of shell effects for superheavy nuclei, the corresponding macro- 
scopic energies as a function of deformation (3 2 with the effective interaction NL3 are pre- 
sented as the dashed lines in Fig. The macroscopic energy is defined as the difference 
between the binding energy and the total shell correction energy which is obtained simi- 
larly as in Section III.C. It can be seen that without the shell effects, the superheavy nuclei 
hardly exist. It is the shell effects that play an essential role to stabilize the superheavy 
nuclei against the fission. In Fig. El the shell correction energies for neutrons, protons, 
and nucleons are respectively presented as the dot-dashed, dashed lines, and solid lines. 
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With the same proton number Z=120 , the shell correction energies E a h e ii in these nuclei 
sensitively depends on the deformation fa- The shell correction energies for protons E p shell at 
I3 2 = are respective -5.5, -4.0, -1.6, -0.7, -1.1, and -1.8 MeV for 292 1 20 , 304 1 20 , 318 1 20, 
348 1 20 , 358 1 20, and 378 1 20, which partly explain the spherical configuration in Fig. [H3 In 
contrast to E P hell , the shell correction energies for neutrons E™ heU at /3 2 = are respective 
—6.1, —2.6, 2.1, —0.4, —5.5, and —13.5 MeV. Due to the large E p shell and / or E™ hell energies, 
the outstanding spherical local minima for 292 1 20 , 304 1 20 and 378 120 are seen in Fig. 
From the constrained RMF calculation, we can conclude that for the Z = 120 isotope chain, 
292 1 20 and 378 1 20 would be doubly magic nuclei with spherical configuration. 

Inspiring by the spherical absolute minimum of 292 1 20 calculated by constrained RMF 
theory with NL3, we use the simple cold fusion model with the doubly magic nuclei 208 Pb 
as a target to estimate the formation cross section a for 292 120 and 304 1 20, although it is 
noticed that the deformed actinide nuclei are recommended to synthesize superheavy nuclei 
with A > 112 The properties of the nuclei including the microscopic energy extracted 
from the constrained RMF calculation are applied in the phenomenological model (3] with 
the reference reaction chosen as 208 Pb( 48 Ca,ln) 255 No. The free parameter C in the model 
is determined by assuming the formation cross section a of the reference reaction, i.e., 
cr(Ci)=500nb and <r(C2)=260nb. In Table HVl the Q value, the height of the static fission 
barrier £?y tat (ER), the optimal excitation energy E* pt , and the formation cross section a at 
E* pt with the input from the constrained RMF for effective interactions NL1, NL3, NLSH, 
and TM1 are listed. The strong dependence of the cross section a on the RMF effective 
interactions can be seen. The largest cross section a is given by NL1. Considered that 
the experimental measure limit of cross section a is above lpb, the two possible doubly 
magic nuclei 292 1 20 and 304 1 20 can be synthesized in the laboratory with the view for NL1. 
However the results from NL1 are not supported by the other RMF interactions. It should 
be mentioned that the results of NL3 are a little different from those of NLSH and TM1, it 
gives nearly accessible cross section (o ~ 1 pb) for 292 120. 

V. SUMMARY 

We have investigated the possible doubly magic nuclei in the superheavy region with 
proton number Z = 100 - 140 and neutron number N= (Z+30) - (2Z+32) within the RCHB 
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theory with interactions NL1, NL3, NLSH, TM1, TW99, DD-ME1, PK1, and PK1R. Shell 
closure are quantified in terms of the two-nucleon separation energies S 2 p and S 2n , the two- 
nucleon gaps S 2p and S 2n , the shell correction energies E p shell and E™ hell , the pairing energies 
Epair an d Epair an d the pairing gaps A p and A n obtained in the RCHB theory. The ex- 
decay half-lives are also discussed. Z=120, 132, and 138 and iV=172, 184, 198, 228, 238, 
and 258 are inferred to be magic numbers. In addition, the spherical configuration of the 
doubly magic nuclei 292 1 20 is supported by examining the potential energy surfaces in the 
deformation-constrained RMF theory and is consistent with other studies applying relativis- 
ts forcesj^J. The shell effects stabilizing the superheavy nuclei are emphasized by extracting 
the shell correction energies form the deformation-constrained RMF calculation. Finally the 
formation cross sections of 292 1 20 and 304 1 20 with 208 Pb as a target are estimated. With 
the view of effective interaction NL1, 172I2O may be synthesized in the current experimental 
setup. 

Acknowledgments 

This work is partly supported by the Major State Basic Research Development Program 
Under Contract Number G2000077407, the National Natural Science Foundation of China 
under Grant Nos. 10025522, and 10221003, and the Doctoral Program Foundation from the 
Ministry of Education in China. 



[1] W. D. Myers and W. J. Swiatecki, Nucl. Phys. A81, 1 (1966). 
[2] H. Meldner, Ark. Fys. 36, 593 (1967). 

[3] S. G. Nilsson, C. F. Tsang, A. Sobiczewski, Z. Szymariski, S. Wycech, C. Gustafson, I.-L. 

Lamm, P. Moller, and B. Nilsson, Nucl. Phys. A131 , 1 (1969). 
[4] U. Mosel and W. Greiner, Z. Phys. 222, 261 (1969). 

[5] P. Moller, J. R. Nix, W. D. Myers, and W. J. Swiatecki, At. Data Nucl. Data Tables 59, 185 
(1995). 

[6] S. Cwiok, V. V. Pashkevich, J. Dudek, and W. Nazarewicz, Nucl. Phys. A410, 254 (1983). 
[7] A. Sobiczewski, Z. Patyk, and S. Cwiok, Phys. Lett. B 224 1 (1989). 



18 



[8] Z. Patyk and A. Sobiczewski, Nucl. Phys. A533, 132 (1991). 

[9] S. Hofmann, et al, Z. Phys. A 350, 277 (1995). 

[10] S. Hofmann, et al, Z. Phys. A 350, 281 (1995). 

[11] S. Hofmann, et al, Z. Phys. A 354, 229 (1996). 

[12] A. Ghiorso, et al, Phys. Rev. C 51, R2293 (1995). 

[13] Yu. A. Lazarev, et al, Phys. Rev. Lett. 73, 624 (1994). 

[14] Yu. A. Lazarev, et al, Phys. Rev. Lett. 75, 1903 (1995). 

[15] Yu. A. Lazarev, et al, Phys. Rev. C 54, 620 (1996). 

[16] Yu. Ts. Oganessian, et al, Phys. Rev. Lett. 83, 3154 (1999). 

[17] Yu. Ts. Oganessian, et al, Phys. Rev. C 62, 041604(R)(2000). 

[18] Yu. Ts. Oganessian, et al, Phys. Rev. C 63, 01 1301 (R) (2000). 

[19] Yu. Ts. Oganessian, et al, Phys. Rev. C 69, 021601(R)(2004). 

[20] P. Moller and J. R. Nix, J. Phys. G 20, 1681 (1994). 

[21] A. Sobiczewski, Sov. J. Part. Nucl. 25, 119 (1994). 

[22] C. L. Wu, M. Guidry, and D. H. Feng, Phys. Lett. B 387,449 (1996). 

[23] G. A. Lalazissis, M. M. Sharma, P. Ring, and Y. K. Gambhir, Nucl. Phys. A608, 202 (1996). 

[24] S. Cwiok, J. Dobaczewski, P.-H. Heenen, P. Magierski, and W. Nazarewicz, Nucl. Phys. A611, 
211 (1996). 

[25] K. Rutz, M. Bender, T. Biirvenich, T. Schilling, P.-G. Reinhard, J. A. Maruhn, and W. 

Greiner, Phys. Rev. C 56, 238 (1997). 

[26] M. Bender, K. Rutz, P.-G. Reinhard, J. A. Maruhn, and W. Greiner, Phys. Rev. C 60, 034304 

(1999). 

[27] J. Meng and N. Takigawa, Phys. Rev. C 61, 064319 (2000). 

[28] W. Long, J. Meng, and S.-G. Zhou, Phys. Rev. C 65, 047306 (2002). 

[29] L. S. Geng, H. Toki, and J. Meng, Phys. Rev. C 68, 061303(R)(2003). 

[30] P. Ring, Prog. Part. Nucl. Phys. 37 193 (1996) 

[31] J. Meng, Nucl. Phys. A635, 3 (1998). 

[32] J. Meng and P. Ring, Phys. Rev. Lett. 77, 3963 (1996). 

[33] J. Meng and P. Ring, Phys. Rev. Lett. 80, 460 (1998). 

[34] J. Meng, I. Tanihata, and S. Yamaji, Phys. Lett. B 419 1 (1998). 

[35] J. Meng, S.-G. Zhou, and I. Tanihata, Phys. Lett. B 532, 209 (2002) 

19 



J. Meng, H. Toki, J. Y. Zeng, S. Q. Zhang, and S.-G. Zhou, Phys. Rev. C 65, 041302(R)(2002). 
R. Smolahczuk, Phys. Rev. C 56, 812 (1997). 

Suk-Joon Lee, J. Fink, A. B. Balantekin, M. R. Strayer, A. S. Umar, P. -G. Reinhard, J. A. 
Maruhn, and W. Greiner, Phys. Rev. Lett. 57, 2916 (1986). 

G. A. Lalazissis, J. Kdnig, and P. Ring., Phys. Rev. C 55, 540 (1997). 
M. M. Sharma and P. Ring, Phys. Rev. C 45, 2514 (1992). 

Y. Sugahara and H. Toki, Nucl. Phys. A579, 557 (1994). 
S. Typel and H. H. Wolter., Nucl. Phys. A656, 331 (1999). 

T. Niksic, D. Vretenar, P. Finelli, and P. Ring., Phys. Rev. C, 66, 024306 (2002). 

W. Long, J. Meng, N. V. Giai, and S.-G. Zhou, |nud-th/03lT03l"| to be published in Phys. 

Rev. C. 

R. Smolahczuk, Phys. Rev. C 59, 2634 (1999). 

H. Kucharek and P. Ring, Z. Phys. A 339, 23 (1991). 

P. Ring and P. Schuck, The Nuclear Many-body Problem 269 (Springer- Verlag, New York, 
1980) 

A. Bohr and B. R. Mottelson, Nuclear Structure Vol. 1 (W. A. Benjamin Inc. , New York, 
Amesterdam, 1969). 

P. Marmier and E. Sheldon, Physics of Nuclei and Particles vol I. (Academic, New York and 
London, 1971). 

Editorial Board of Nucl. Phys. Rev., Nucl. Phys. Rev. vol. Nuclear Chart (Science Press, 
Lanzhou, 2003). 

M. Bender , W. Nazarewicz, and P.-G. Reinhard, Phys. Lett. B. 515, 42 (2001). 
V. E. Viola Jr. and G. T. Seaborg, J. Inorg. Nucl. Chem. 28 , 741 (1966). 
P. Armbruster, C. R. Physique 4, 571 (2003). 



20 



TABLE I: The effective interactions PK1 and PK1R 44 1. The mass of the isovector vector meson 
is taken as m p =763 MeV for PK1 and PK1R. 



Effective 


m a m u g a g w 


9 P 


92 


53 c 3 d 3 


interactions 


[MeV] [MeV] 


[fm- 1 ] 






PK1 


514.089 784.254 10.322 13.013 


4.530 


-8.169 


-9.998 55.636 


PK1R 


514.087 784.222 10.322 13.013 


4.550 


-8.156 


-10.198 54.446 350 
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TABLE II: The possible magic proton number suggested by two-proton separation energies, two- 
proton gaps, shell correction energies, pairing energies and effective pairing gaps for proton with 
interactions NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 



Effective 
interactions 


Quantity 


Magic proton number candidates 


106 114 120 126 132 


138 


NL1 


&2p 

hp 

shell 
E P ■ 

pOflT 

A p 


■ ■ ■ ■ 

■ ■ ■ ■ ■ 

■ ■ ■ ■ ■ 

■ ■ ■ ■ ■ 




S2p 


■ 1 


1 ■ 


■ 




hp 








NL3 


E p 

shell 


■ 1 


1 ■ 


■ 




E p ■ 


■ ■ 1 


1 ■ 


■ 




pair 










A p 


■ ■ 1 


1 ■ 


■ 




S2p 


■ 1 


1 ■ 


■ 




hp 


■ ■ 1 






NLSH 


^ shell 


■ ■ 1 


1 ■ 


■ 




fpP 


■ ■ 1 








pair 










A p 


■ ■ 1 








S2p 










hp 








TM1 


E p 

shell 










E p ■ 










PCLIT 










Ap 










S2p 










hp 


■ ■ 1 






TW99 


E p 

shell 


■ ■ 1 








E p • 


■ ■ 1 


1 ■ ■ 






pair 










Ap 


■ ■ 1 


1 ■ ■ 






S2p 










hp 








DD-ME1 


E p 

shell 










E p ■ 


■ ■ 1 


1 ■ ■ 






pair 










Ap 


■ ■ 1 


1 ■ ■ 




PK1 


J2p 

hp 

^ shell 

E p ■ 

pair 

Ap 


■ ■ ■ ■ ■ 

■ ■ ■ ■ ■ 

■ ■ ■ ■ ■ 


PK1R 


S2p 

hp 

^shell 

E P ■ 

pair 

A p 


■ ■ ■ ■ ■ 

■ ■ ■ ■ ■ 

■ ■ ■ ■ ■ 
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TABLE III: The possible magic neutron number suggested by two-neutron separation energies, 
two-neutron gaps, shell correction energies, pairing energies and effective pairing gaps for neutron 
with interactions NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 



Effective 


Quantity 


Magic neutron number candidates 


interactions 


138 164 172 184 198 216 228 238 252 258 274 




$2n 


■ ■■■ ■■■■■ 




fori 


■ ■■■ ■■ ■■ 


NL1 


rpn 
shell 


■ ■ 




E n ■ 

pair 

A n 


■ ■■■■ ■■■■■ 

■ ■■■■ ■■■■■ 




£>2n 


■ ■■■ ■■■■■ 




fori 


■ ■■■■ ■■ ■■ 


NL3 


rpn 

^shell 






E n ■ 

pair 

A n 


■ ■■■■■■■■■■ 

■ ■■■■■■■■■■ 




S2n 


■ ■■■■■■■■■■ 




fori 


■ ■■■■■■■ ■■ 


NLSH 


T?n 

^shell 


■ ■ ■ ■ ■ 




rpn 

pair 

A n 


■ ■■■■■■■■■■ 

■ ■■■■■■■■■■ 




&2n 


■ ■■■ ■■■■■ 




fori 


■ ■■■ ■■ ■■ 


TM1 


rpn 

^shell 
rpn 
pair 

A n 


■ ■■■■ ■■■■■ 




&2n 


■ ■■■■■■■■■■ 




fori 


■ ■■■■■■■ ■■ 


TW99 


rpn 

^shell 


■ ■ ■ ■ ■ 




rpn 

pair 


■ ■■■■■■■■■■ 

■ ■■■■■■■■■■ 




SWi 


■ ■■■■ ■■■■■ 




fori 


■ ■■■■ ■■ ■■ 


DD-ME1 


rpn 

^shell 

E n ■ 

pair 

A n 


■ ■ ■ ■ 

I I I I I I I I I I 




02n 






fori 


■ ■■■■■■■ ■■ 


PK1 


rpn 

^shell 






E n ■ 

pair 

A n 


■ ■■■■■■■■■■ 

■ ■■■■■■■■■■ 




S2n 


■ ■■■■■■■■■■ 




fori 


■ ■■■■■■■ ■■ 


PK1R 


rpn 

^shell 






E n ■ 

pair 


■ ■■■■■■■■■■ 

■ ■■■■■■■■■■ 

1 c*=) 



TABLE IV: The Q value, the height of the static fission barrier £?|* a *(ER), the optimal excita- 
tion energy E* pt and C-dependent formation cross section a at E* pt for the cold fusion reaction 
208 Pb( 85 Sr,ln) 292 120 and 208 Pb( 97 Sr,ln) 304 120 with the input from the constrained RMF for effec- 
tive interactions NL1, NL3, NLSH and TM1 



Cold fusion reactions 


RMF 


Q 


Bf at (ER) 


TP* 

^opt 


a 




Interactions [MeV] 


[MeV] 


[MeV] 


aid) 


a(C 2 ) 




NL1 


296.51 


10.62 


14.78 


15.8pb 


9.5pb 


208 Pb( 85 Sr,ln) 292 120 


NL3 
NLSH 


302.84 
297.88 


8.15 
7.21 


14.92 
13.94 


l.lpb 
0.26fb 


720fb 
0.16fb 




TM1 


296.83 


5.74 


12.83 


0.47fb 


0.29fb 




NL1 


308.42 


12.62 


9.96 


86.8nb 


55.1nb 


208 Pb( 97 Sr,ln) 304 120 


NL3 
NLSH 


305.09 
300.37 


3.98 
5.15 


9.55 
10.49 


0.64fb 
13.8fb 


0.40fb 
8.6fb 




TM1 


299.73 


5.51 


11.21 


20.5fb 


12.6fb 




FIG. 1: Nuclear chart in the superheavy region. The long-dashed and solid lines represent the /te- 
stability lines by formula N—Z = 6x 10~ 3 -A 5 / 3 kgf and Z = ^ kifl respectively. 

^ J 1.98 + 0.0155 • A 2 / 3 w 

The dots lying at the lower-left corner represent the superheavy nuclei observed or declared to be 

observed experimentally. 
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FIG. 3: The two-neutron separation energies S2n(N, Z) = Eb(N, Z) — Eb(N — 2, Z) as a function 
of mass number A obtained by RCHB calculation with effective interactions NL1, NL3, NL-SH, 
TM1, TW-99, DD-ME1, PK1, and PK1R, respectively. 
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FIG. 4: The two-proton gaps 52 P (N, Z) = 5 > 2 P (A r , Z) — S2 P (N, Z + 2) as a function of proton number obtained by RCHB calculation with 
effective interactions NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1, and PK1R, respectively. 
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FIG. 5: The two-neutron gaps 62 n (N, Z) = S2n{N, Z) — S2n(N + 2, Z) as a function of neutron number obtained by RCHB calculation with 
effective interactions NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1, and PK1R, respectively. 




100 110 120 130 140 100 110 120 130 140 

Proton Number Z Proton Number Z 



FIG. 6: The shell correction energies for proton E?,,, as a function of proton number obtained by RCHB calculation with effective 
interactions NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 




FIG. 7: The shell correction energies for neutron E™ hell as a function of neutron number obtained by RCHB calculation with effective 
interactions NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 
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FIG. 8: The pairing energies for proton E^ air as a function of proton number obtained by RCHB calculation with effective interactions NL1, 
NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 
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FIG. 9: The effective pairing gaps for proton A p as a function of proton number obtained by RCHB calculation with effective interactions 
NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 
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FIG. 10: The pairing energies for neutron E™ air as a function of neutron number obtained by RCHB calculation with effective interactions 
NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 
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FIG. 11: The effective pairing gaps for neutron A n as a function of neutron number obtained by RCHB calculation with effective interactions 
NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1 and PK1R, respectively. 
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FIG. 12: The a-decay half-lives T a as a function of neutron number obtained by RCHB calcu- 
lation with effective interactions NL1, NL3, NL-SH, TM1, TW-99, DD-ME1, PK1, and PK1R, 
respectively. 
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FIG. 13: The binding energy and the equivalent macroscopic energy of 292 > 304 > 318 > 348 > 358 ' 378 i20 
calculated in the constrained RMF theory with effective interaction NL3. 
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FIG. 14: The shell correction energy of ^92,304,ji8,348,J58,378 12 o calculated in the constrained RMF 
theory with effective interaction NL3. 
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